NASA Technical Memorandum 1006H 


THERMAL STRESS IN HIGH TEMPERATURE 
CYLINDRICAL FASTENERS 

(NASA-TM-10061 1) THERMAL STRESS IN HIGH N88-22438 

TEMPERATURE CYLINDRICAL FASTENERS {NASA) 

13 p CSCL 20K 

Unclas 
G3/39 0141241 


Max L. Blosser 


MAY 1988 


NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665-5225 



THERMAL STRESS IN HIGH TEMPERATURE CYLINDRICAL FASTENERS 


Max L. Blosaer 

NASA Langley Research Center 
Hampton. VA 23665 


Abstract 

Uninsulated structures that are fabricated 
from carbon or silicon-based structural materials 
and that are allowed to become hot during flight 
are attractive for the design of some components 
of hypersonic vehicles. They have the potential 
to reduce weight and Increase vehicle efficiency. 
Because of manufacturing constraints, these 
structures will consist of parts which must be 
fastened together. The thermal expansion mismatch 
between conventional metal fasteners and carbon or 
silicon-based structural materials may make it 
difficult to design a structural Joint which is 
tight over the operational temperature range 
without exceeding allowable stress limits. In 
this study, algebraic, closed-form solutions for 
calculating the thermal stresses resulting from 
radial thermal expansion mismatch around a 
cylindrical fastener are developed. These 
solutions enable a designer to quickly evaluate 
many combinations of materials for the fastener 
and structure. Using the algebraic equations 
developed in this study, material properties and 
Joint geometry were varied to determine their 
effect on thermal stresses. Finite element 
analyses were used to verify that the closed-form 
solutions derived in this study give the correct 
thermal stress distribution around a cylindrical 
fastener and to investigate the effect of some of 
the simplifying assumptions made in developing the 
closed-form solutions for thermal stresses. 


Introduction 

The National Aerospace Plane (NASP) program 
has generated a renewed interest in research 
relating to hypersonic vehicles. The severe 
aerodynamic heating encountered by reusable 
hypersonic vehicles makes the structural design of 
such vehicles a challenge. Uninsulated 
structures, which are alloweq to become hot 
during flight, are attractive for the design of 
some components of hypersonic vehicles because 
they have the potential to reduce weight and 
increase vehicle efficiency. Some hot structures 
being considered may reach temperatures 
approaching 3000 # F. Carbon or silicon-based 
materials will likely be useful for such 
structures because they are light and retain 
strength at elevated temperature. 

An example of such a hot structure is a 
carbon-carbon body flap which was studied as a 
replacement for the insulated aluminum flap on the 
Shuttle (ref. 1). In addition, NASA Langley 
Research Center has a current contract 
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("Development of Oxidation Protected Carbon-Carbon 
Lightly Loaded Structures" - contract NASI -1 8462) 
which will include the design and fabrication of a 
carbon-carbon control surface test article for 
NASP, such as that shown in figure 1. Because of 
manufacturing constraints, these structures and 
similar hot structures will consist of smaller 
pieces which must be fastened together. Ductile 
metal fasteners are desirable for assembling these 
multiple pieces. However, the thermal expansion 
mismatch between the metal fastener and carbon or 
silicon-based structural materials makes it 
difficult to design a structural Joint which is 
tight over the operational temperature range 
without exceeding allowable stress limits. 



Fig. 1 : Hot structures for reusable hypersonic 
vehicles 


The problems encountered using a standard 
cylindrical fastener in a high temperature Joint 
are illustrated in figure 2. In the figure the 
fastener is assumed to have a much larger 
coefficient of thermal expansion (CTE) than the 
structure. If* the fastener is snug in the radial 
direction at room temperature, high thermal 
stresses will result at elevated temperature. If 
sufficient radial clearance is left around the 
fastener so that it is Just snug at elevated 
temperature, then it may be too loose at room 
temperature. Similarly, if the fastener is 
axially snug at room temperature it will be loose 
at elevated temperature, and it must be 
prestressed considerably at room temperature to 
remain snug at elevated temperature. The radial 
thermal stress or axial, prestress required to 
maintain an acceptably tight Joint over the 
operational temperature range may cause premature 
failure of the Joint. Therefore, a Joint 


consisting of a standard cylindrical fastener that 
Joins pieces of the selected high temperature 
material should be analyzed to determine if there 
are satisfactory combinations of radial clearance 
and axial prestress at room temperature which will 
produce acceptable stresses at elevated 
temperature. 




1 Joint condition 1 


Room temporatura 

High temperature 

Racial (fraction 

Snug 

Loose 

Hgh thermal ttsesee 
Snug 

Axial direction 

Snug 

Pre- stressed 

Loose 

Snug 


Fig. 2: Problem areas for cylindrical fastener 
Joints in high temperature applications 


In this paper, simple solutions for 
calculating the thermal stresses resulting from 
radial thermal expansion mismatch around a 
cylindrical fastener are developed to enable a 
designer to evaluate quickly many combinations of 
materials. The solutions developed do not address 
the axial prestress. 

Although similar solutions were found in the 
literature, no explicit solution for the thermal 
stresses around a cylindrical fastener was found. 
References 2— *» contain axisymmetric plane-stress 
solutions, but these solutions do not involve 
temperature or thermal expansion. References 5-9 
contain similar solutions involving thermal 
expansion, but do not address quite the same 
problem as this paper. ( A more detailed 
discussion of references 2-9 is given in Appendix 
A.) Consequently, algebraic equations for the 
thermal stresses in a finite ring of material 
around a cylindrical fastener of a different 
material are developed in this study. Similarly, 
algebraic equations for the thermal stresses 
around a hollow cylindrical fastener are also 
developed. Equations are also presented that 
account for the effect of an initial clearance 
between the fastener and surrounding material. 

The equations developed in this study are used to 
study the effect of varying Joint geometry and 
material properties on the thermal stresses in a 
cylindrical fastener Joint. Three-dimensional 
finite element analyses are used to verify these 
solutions and to evaluate the effect of some of 
the simplifying assumptions on the magnitude and 
distribution of the thermal stress. 


Symbols 

a radius of fastener 

a t inside radius of hollow fastener 

A,B,C,D constants defined by boundary conditions 
b radius of structure around fastener 

E modulus of elasticity 

L length 

L # length at initial temperature 

n integer 

P pressure at interface between fastener 

and structure 

r,6,z cylindrical coordinates 
T temperature 

T, initial temperature 

u displacement 

a coefficient of linear thermal expansion 

(CTE) 

5 initial clearance around a cylindrical 

fastener 

+ Airy stress function 

c strain 

v Poisson's ratio 

a stress 

Subscripts 

f fastener 

8 Sap 

max maximum 

r r direction 

s structure around the fastener 

z z direction 

8 8 direction 


Closed-Form Solutions 


The objective of the following analysis is to 
develop a simple equation to predict the thermal 
stresses around a cylindrical fastener at elevated 
temperature. The thermal stresses predicted by 
this analysis occur because the fastener has 
greater thermal expansion in the radial direction 
than the surrounding structure. The fastener is 
also assumed to have greater thermal expansion in 
the axial direction than the surrounding 
structure and, with increasing temperature, will 
tend to loosen axially, relieving any axial 
prestress. If appreciable normal stresses develop 
at the interface between the fastener and 
structure, friction will inhibit the axial 
slippage of the fastener and thereby induce out- 
of-plane stresses in the Joint. However, in the 
present analysis the assumption is made that there 
is no friction between the fastener and 
surrounding structure. Thus, no significant 
stresses are introduced in the thickness direction 
,of the material and a condition of plane stress 
results. The fastener and surrounding structure 
are assumed to undergo a uniform temperature 
change and to expand radially as isotropic 
materials although the CTE's of the two materials 
differ. For simplicity-, the fastener is assumed 
to be surrounded by a disk of material with the 
fastener at the center. In most Joints the shape 
of the structure will be much more complicated, 
and the thermal stresses will vary somewhat around 
the fastener. Also, for simplicity the ends of 
the bolt, including the head, are neglected (see 
fig. 3). With these assumptions the problem 
becomes an axisymmetric, plane-stress problem. 
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Fig. 3: Idealization of joint geometry 


Solid Fastener Solution 


The general equations developed in Appendix 
B (B12 and B13) for axlsymmetric stresses and 
displacements are now used to find the thermal 
stresses around an idealized cylindrical fastener 
(fig. 3) at elevated temperature. If the fastener 
and structure are assumed to be at a uniform 
temperature and to expand relative to the center 
of the fastener, the thermal displacement is 
simply proportional to the radial position. 
Therefore, by superposition a thermal expansion 
term la added to equation (B13) for radial 
displacement to produce: 


For the structure surrounding the fastener, the 
boundary conditions are: 

at r ■ a, a (a) - -P and at r - b, o (b) - 0 
r a s 


Applying these boundary conditions to equations 
(B12), the constants A and C are found to be: 


A - 


P b 


bi2 


-(f) 


and C - - % ( 1 v- = ) (5) 

2 1 - (f) 2 


The expression for radial displacement in the 
structure surrounding the fastener therefore 
becomes: 


-I* 


» 1 - (f) 


w-oHr Cl ) ♦ (1-v )r) ♦ r a AT 

r 9 9 9 


(6) 


Displacement compatibility at the Interface 
between the fastener and structure, r - a, 
requires that the radial displacements be equal at 
that point. Equations (3) and (6), with r • a, 
can be combined and solved for the interface 
pressure, P. 



E a ((f) 2 -1)( « f - «, ) AT 
E„ h 5 

(Uv 8 ) ♦ (1-v a ) ♦ ^ ((f) -1) (1-V f > 


(7) 


The stress distributions in the structure 
surrounding the fastener can be found by combining 
equations (B12) and (5) to produce: 


u r “ f ( - ♦ 2 C (1-v) r ) ♦ r a AT (1) 

The stresses and displacement of the fastener 
and surrounding material are determined by 
applying the appropriate boundary conditions. For 
the fastener: 

at r - a, a (a) - -P 
r r 

where P is the pressure at the Interface between 
the fastener and structure and 


< s > 2 - 1 

V rl #TT 


A 2 . i 

V p 1 <j> 2 - 1 


where P is given by equation (7). 


(8) 


(9) 


at r - 0, o_ (0) is finite. 
r f 

Applying these boundary conditions to equations 
(B12), the constants A and C are found to be: 

A - 0, and C - - f (2) 


Therefore, the radial displacement for the 
fastener is given by: 


P r 


(1-v f ) 


♦ r a. AT 


(3) 


Substituting equations (2) into (B12) shows that 
the fastener is in hydrostatic compression, 



Several limiting cases of the solution for 
thermal stresses around a cylindrical fastener are 
now considered. These limiting cases give insight 
into the behavior of the Joint and provide simpler 
expressions for thermal stresses in the Joint. If 
the structure around the fastener is shrunk to an 
infinitely thin ring, the stresses in the fastener 
approach zero. The only nonzero stress is the 
hoop stress in the ring of structure which is 
given by the simple expression: 

°8 ’ E a (a f ' V AT (10) 


If the fastener is surrounded by an infinite 
sheet (b/a ♦ •), equation (7) reduces to the 
following form: 


P - 


E a (tt f - a a ) AT 


(1 ♦ v fl ) 



V 


(11) 
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and the stresses in the structure around the 
fastener become: 


P. 


a 



( 12 ) 


If the fastener material Is much stlffer than 
the surrounding structure ( E f » E a ), then 

equation (11) reduces to: 


P 


E 

s 


(®f - a a ) AT 


1 
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V 


a 


(13) 


E s ((^) 2 - i)(q f -c 3 ) AT 


(J) 2 (Uv a ) + (1-v a )* J 

f a^- a 

(17) 


The stresses in the structure surrounding the 
fastener can be found /by substituting the value of 
P from equation (17) into equations (8) and (9). 
The stress distributions in the hollow fastener 
can be found by combining equations (B12), and 
(15) to produce: 


and the stresses may be found by substituting 
equation (13) into equations (8) and (9). 
Similarly, if the surrounding structure is much 
stiffer than the fastener material ( E g >> E f ), 

then equation (11) reduces to: 


- P ( 


(=i) 2 - 1 


< r > 2 - 1 


P 


E r {a f " a s ) AT 

i - V. 


Hollow Fastener Solution 


(HO 


o a ■ P ( 


(-i) 2 ♦ 1 


< r > 2 - 1 


(18) 


(19) 


The solution for thermal stresses of a hollow 
cylindrical fastener can be found by following a 
similar procedure. The equations for stresses and 
displacements in the structure surrounding the 
fastener will have the same form for the hollow 
fastener as for the solid fastener, however, the 
equation for the interface pressure between' the 
fastener and surrounding structure will be 
different. 

The stresses and displacement for the 
idealized hollow fastener (fig. 3) can be found by 
applying the following boundary conditions: 

at r • a, , o_(a. ) - 0 

1 r l 

and at r - a , <r (a) - -P 
r 


Fastener With Initial Clearance 

Fasteners may have an initial clearance at 
room temperature. Equation (7) may be easily 
modified to account for an initial clearance. The 
temperature change which closes the gap due to the 
initial clearance must be calculated. This 
temperature can be found by equating the radial 
position of the expanded fastener and hole. As 
discussed in reference 10, if the CTE's are 
assumed to be independent of temperature, the 
exact expression for thermal expansion can be 
found by integrating the equation: 



dL - L. o dT 

(20) 

to produce 

L - L. .“ <T - T -> 

(21) 


Applying these boundary conditions to equations 
(B12), the constants A and C are found to be: 


A - 


„ 2 2 
P a a^ 

~2 2 

* - \ 


and 


P a 


2 (a 2 - a 2 t ) 


(15) 


Equating the radial position of the expanded 
fastener and hole therefore produces: 


l «r AT, ) 


A 6 


( 

(a ♦ «) e 


“s AT g 


( 22 ) 


Solving equation (22) for the temperature change 
produces : 


The expression for radial displacement in the 
fastener, therefore, becomes: 


£-4 )lir-(l+v f ) + (1-v f ) r) ♦ r « f AT 

f 1 - (— f 

v a ’ 

(16) 


As in the solid fastener solution, 
displacement compatibility at the interface 
between the fastener and surrounding structure, 
r ■ a, requires that the radial displacements be 
equal at that point. Equations (6) and 06), with 
r • a, can be combined and solved for the 
interface pressure P to give: 


AT. 


ln l' * I > 

“f • “s 


(23) 


where 6 is the initial clearance. The 
temperature change calculated using equation (23) 
can be subtracted from the total temperature 
change, and the net temperature ohange can be 
substituted into equation (7) to find the 
corresponding thermal stresses. Note, however, 
that the term (AT - ATg) must be positive. 


Finite Element Analyses 


The finite element method is commonly used to 
numerically analyze the stresses in Joints. 
Therefore, finite element solutions were obtained 
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for comparison with the mathematically derived 
solutions presented In this study. Finite element 
analyses were used to verify that the equations 
derived In this study give the oorrect thermal 
stress distribution around a cylindrical fastener. 
In addition the finite element analyses were used 
to investigate the effect of friction and the 
protruding end of the fastener on the thermal 
stress distribution. An existing, general-purpose 
finite element program. Engineering Analysis 
Language (EAL) (ref. 11), was used for the 
analyses. 

A three-dimensional finite element model 
(see figure 4) was used to represent a solid 
cylindrical fastener Joint configuration. A five- 
degree wedge-shaped model of the fastener and 
structure was used as indicated on the figure. 
Similar models were used to model both a hollow 
fastener Joint and a Joint with a fastener which 
protruded from the structure. 



v - 4 

« - 1 xW* kw/kv/V 


Fig. it: Finite element model of solid cylindrical 
fastener and surrounding structure 

Particular, but arbitrary, dimensions were 
chosen for each Joint modeled. The radius of the 
cylindrical fastener was 0.1 Inch and the radius 
of the surrounding structure was 1 .0 Inch. The 
thickness of the cylindrical fastener model was 
0.15 Inch, but the symmetry constraints applied on 
the lower surface made the model respond as though 
it were 0.3 Inch thick. The Inside radius of the 
hollow fastener was 0.05 inch. For the model of 
the protruding fastener, the fastener extended 0.1 
inch above the structure. An arbitrary set of 
Isotropic stiffness properties was used (see 
figure A). 

Three-dimensional, six-node and eight-node, 
linear elements were used to model the fastener 
and eight-node linear elements were used to model 
the surrounding structure. The radial dimensions 
of the elements were varied parabollcally to 
achieve a finer mesh near the Interface between 
the fastener and surrounding structure. The finer 
mesh near the Interface was required for efficient 
prediction of the expected stress gradients In 
that region. The resulting finite element mesh 
can be seen In figure M. Nodes on the sides and 
bottom of the model were constrained to remain In 
their respective planes but were free to move 
within those planes. 


The most challenging portion of the finite 
element analysis was modeling the contact between 
the fastener and surrounding material along the 
bearing surface. Adjacent nodes along this 
boundary were connected by zero-length elements 
which had high stiffness perpendicular to the 
bearing surface and no stiffness tangent to the 
surface. The elements allowed relative motion 
along the boundary between the adjacent nodes, but 
not perpendicular to It. The model was set up so 
that If one of these elements was found to be in 
tension, the element stiffness could be reduced to 
an Insignificant value. For the cases analyzed, 
however, none of the zero length elements was in 
tension. For the solid fastener, one case was 
analyzed with the zero length elements having 
large (essentially infinite) stiffnesses In all 
degrees of freedom. This case simulated the 
highest possible stresses due to friction. 


Results and Discussion 

Comparison of Closed-Form and Finite Element 
Solutions 

The thermal stresses calculated from the 
olosed-form solutions developed In this paper are 
oompared to those calculated from the finite 
element models previously described. Figure 5 
shows the thermal stress comparison for a solid 
fastener subjected to a 1000 *F temperature 
Increase. The radial and hoop stresses are shown 
as a function of radial position for both the 
fastener and structure. The shaded region of the 
figure represents the fastener, and the unshaded 
region represents the surrounding structure. 
Stresses calculated using the three-dimensional 
finite element model are Indicated by symbols, and 
the plane-stress solution stresses are Indicated 
by solid lines. The results show excellent 
agreement for both the magnitude and distribution 
of thermal stresses. The peak stress of 
approximately 80,000 psl Is Indicative of the 
large thermal stresses which may occur In hot 
structures with cylindrical fasteners. The peak 
stresses In the structure occur at the Interface 
between the fastener and structure. The magnitude 
of the oompressive radial stress Is slightly less 
than the magnitude of the tensile hoop stress. 



Fig. 5: Thermal stresses for solid cylindrical 
fastener 
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The fastener Is in hydrostatic compression with 
stress equal to the peak radial stress in the 
structure, as shown by the constant hoop and 
radial stresses in the fastener. Because the 
struct ural ho op and radial stresses are tensile 
and compressive respectively, the von Hises 
stress in the structure will be considerably 
larger than in the fastener. The geometry and 
loading are axisymmetric, so the in-plane shear 
stress is Inherently zero. The plane-stress 
solution neglects the z stress and out-of-plane 
shear stresses. In the absence of friction 
between the fastener and structure, the three- 
dimensional finite element analysis also predicts 
negligible out-of-plane stresses for the 
relatively thick model analyzed. 

Figure 6 shows a similar comparison of 
thermal stresses for a hollow fastener. The 
geometry, material properties, and loading are the 
same as for figure 5, except that the fastener is 
hollow. The shaded band represents the wall of 
the hollow fastener. Again, the only significant 
stresses are the radial and hoop stresses in the 
fastener and structure. The figure shows the 
excellent agreement between the three-dimensional 
finite-element and plane-stress solutions. The 
stress distribution in the structure is similar to 
that for the solid fastener, however, the 
magnitudes of the stresses are lower. The stress 
distribution in the hollow fastener is completely 
different than in the solid fastener. The radial 
stress goes to zero on the inner surface of the 
fastener, as expected. The hoop stress in the 
hollow fastener is much larger than that in the 
solid fastener. The peak hoop stress occurs on 
the inner surface of the hollow fastener, and for 
this example, is more than double the stress in a 
solid fastener. Again, the out-of-plane stresses 
predicted by the three-dimensional finite element 
solution were negligible. 



Fig. 6: Thermal stresses for hollow fastener 


Effects of Geometry 
Thermal Stress 


and Material Properties on 


The equations for thermal stress around a 
solid cylindrical fastener were used to determine 


the effect of varying material properties and 
geometry. Figure 7 shows the maximum thermal 
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Fig. 7: Maximum structural thermal stresses for 
solid fastener as a function of ratio of 
structure-to-fastener radius 


stresses in the structure around a solid fastener 
as a function of the ratio of the structure radius 
a to the fastener radius b . The thermal 
stress is nondlmenslonallzed by dividing by an 
expression that represents the maximum possible 
thermal stress in the structure. The ratio of 
structural modulus to fastener modulus la varied 
to produce a set of curves on the figure. The 
ratio b/a was varied between one and ten. For b/a 
> 5, the maximum stresses become almost 
Ins ensitive to further changes in b/a an d the 
simpler equations (11) and (12) can be used to 
calculate thermal stresses. As b/a approaches 
unity, the structure becomes a thin ring around 
the fastener and the radial stresses approach 
zero. The ring of material cannot restrain the 
greater thermal expansion of the fastener, and 
almost all of the thermal expansion mismatch is 
accommodated by hoop strain in the ring. The non- 
dimensional hoop stress therefore approaches unity 
(see eqn. (10)). When the modulus of the 
structure is greater than that of the fastener, 
more of the thermal expansion mismatch is 
accommodated by strain in the fastener, and the 
nondimens ional stress in the structure is reduced. 
Conversely, if the modulus of the structure is 
less than that of the fastener, the nondi mens ional 
stress in the structure is increased. 

The Poisson's ratios of the fastener and 
structure were also varied to determine their 
effect on thermal stress. In figure 8 the maximum 
thermal stress (hoop stress at r-a) in a solid 
fastener in an infinite sheet of structure is 
shown as a function of Poisson's ratio of the 
fastener. The Poisson's ratio of the structure is 
also varied to produce a set of curves. The 
values of Poisson's ratio are varied from zero to 
0.5, the range possible with isotropic materials. 
Increasing the fastener 'Poisson's ratio 
effectively stiffens the fastener, which is in 
two-dimensional hydrostatic compression, and 
thereby Increases the peak thermal stress. 
Increasing the structure Poisson's ratio, however, 
decreases peak thermal stresses. 
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Fig. 8: The effect of Poisson's ratio on solid 
fastener maximum thermal stress 


Figure 9 shows a comparison of the maximum 
thermal stresses of a hollow fastener to those of 
a solid fastener with similar geometry, material 
properties, and loading. The ratio of hollow-to- 
solid fastener stresses Is shown as a function of 
the ratio of the hollow fastener lnner-to- outer 
radius. The figure shows that even a small hole 
in the center of a fastener increases the stress 
In the fastener by a factor of 2. The stress 
increase occurs because the stress distribution Is 
changed by the Introduction of the hole. On the 
inner surface of the fastener the radial stress 
must go to zero so that the fastener cannot be in 
hydrostatic compression. The hoop stress 
increases significantly, compared to a solid 
fastener. As the wall of the hollow fastener 
becomes thinner, the stress In the structure is 
reduced, but the fastener stress remains high. 

When the fastener modulus is greater than that of 
_the structure, the fastener stress ratio Increases 
as the fastener wall is made thinner. When the' 
moduli are equal, the fastener stress ratio 
remains constant with wall thickness. When the 
structural modulus is greater, the fastener stress 
ratio decreases as the fastener wall Is made 
thinner. For a specified wall thickness, a 
stiffer fastener modulus produces a lower stress 
ratio. 



Fig. 9: Thermal stresses In hollow fastener 
compared to those in a solid fastener 


Fastener Protrusion Beyond Surrounding Structure 

Host fasteners extend through the structure 
being fastened. However, to develop the algebraic 
equations in this study, the fastener was assumed 
not to protrude beyond the surrounding structure. 

A finite element model was used to investigate the 
effects of this assumption by calculating the 
thermal stresses for a particular fastener which 
protrudes beyond the structure. Figure 10 shows 
the distribution of radial stress In such a joint. 
Away from the fastener the stresses are uniform 
through the thickness, which is consistent with 
the plane-stress solution. Near the fastener the 
stresses In the structure increase slightly toward 
the protruding end of the fastener, compared to 
the stresses in a Joint in which the fastener does 
not protrude. The portion of the fastener 
protruding beyond the structure is free to 
thermally expand radially, but the portion of the 
fastener within the structure is constrained by 
the surrounding structure. Thus, the fastener 
expands sore near the protruding end of the 
fastener than towards the center of the structure. 
The presence of a fastener head would therefore be 
expected to Increase the local, peak thermal 
stresses in the structure even more than a 
protruding fastener shank. The peak stresses for 
a Joint with a protruding fastener are shown In 
Table 1 , along with the stresses predicted using 
the plane-stress solution developed in this study. 
The peak radial stress In the structure is about 
15 percent higher for the protruding fastener. 
However, these higher stresses are localized at 
the fastener/structure Interface at the free 
surface and for many materials may be alleviated 
by local yielding. Significant out-of-plane 
stresses, not present in the plane-stress 
solution, were calculated for the Joint with a 
protruding fastener. However, these out-of -plane 
stresses are low compared to the in-plane 
stresses. If brittle materials or materials with 
low out-of-plane strengths are used In the Joint, 
a more detailed analysis of the Joint may be 
required. 



Fig. 10: Radial thermal stress distribution (psl) 
for solid fastener protruding beyond 
structure 
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Table 1: Maximum calculated thermal stresses (pal) Fig . 11s Temperature change required to close 

initial clearance 


Friction Between Fastener and Structure 

Finite element analysis was also used to 
investigate the maximum possible thermal stresses 
due to friction. The adjacent nodes on the 
boundary between the structure and fastener were 
forced to remain together. The maximum stresses 
calculated for the Joint are shown in Table 1 . 

The peak stresses for the worst case friction are 
almost double the stresses calculated assuming no 
friction. The largest shear stresses occur along 
the interface between the fastener and structure. 
The shear stresses Increase linearly along the 
length of the interface and reach a maximum near 
the free surface. In an actual Joint some 
slipping will likely occur along the Interface to 
alleviate some of the frictional stresses. This 
finite element analysis and the plane-stress 
solution represent the upper and lower bounds, 
respectively, of the thermal stresses for this 
particular idealized Joint. The large difference 
in stresses indicates that friction can have a 
significant effect. Therefore, further 
investigation may be necessary to determine the 
effect of friction on the thermal stresses around 
a cylindrical fastener. A designer might also 
consider taking steps to reduce friction in a 
Joint, such as using a high temperature lubricant. 

Fastener With an Initial Clearance 

The temperature change required to close an 
Initial clearance around a cylindrical fastener is 
shown as a function of clearance to diameter ratio 
in figure 11. A family of curves is shown for 
various differences in thermal expansion 
coefficients between fastener and structure. As 
expected, for larger differences in thermal 
expansion coefficients, smaller increases in 
temperature will close the initial clearance. 

Table 2, from reference 12, lists the clearances 
for several categories of fits. On figure 11, the 
ranges for close-sliding fit, medium running fit, 
and loose running fit are indicated for fasteners 
with 0.2 and 1 .0 inch diameters. For a given type 
of fit, a smaller diameter fastener will require a 
larger temperature increase to close the initial 
clearance. Consequently, for a given type of 
sliding or running 'it, a larger fastener will 
tend to have higher thermal stresses than a 
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Table 2: Limits of clearance for running and 
sliding fits from ref. 12 


smaller fastener. The effect of an initial 
clearance can be seen by considering the sample 
Joint analyzed previously In this study, with 
geometry and material properties defined in figure 
M. For a 1000 *F temperature Increase, the 
thermal stresses around a fastener with a close- 
sliding fit will be reduced by between eight and 
28 percent compared to a fastener with no initial 
clearance. For a medium fit the thermal stress 
reduction would be reduced by between 44 and 100 
percent, and for a loose fit there would be no 
thermal stress. For the same conditions with one 
inch diameter fastener, the thermal stresses for 
the close fit would be reduced by between three 
and 11 percent; the stresses for the medium fit 
would reduced by between 18 and 36 percent; and 
the stresses for the loose fit would be reduced by 
between 50 and 100 percent. 


Concluding Remarks 

Algebraic equations derived from closed-form 
solutions of elasticity equations are presented 
for thermal stresses in solid and hollow 
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cylindrical fasteners. The fastener Mas idealized 
as a pin surrounded by an annulus of structural 
material so that an ax i symmetric plane-stress 
solution could be developed. Friction between the 
fastener and structure was neglected. For these 
assumptions, excellent agreement was obtained 
between the thermal stresses predicted by the 
algebraic equations and by a three-dimensional 
finite element analysis. 

An equation was also developed to predict the 
temperature required to close an initial clearance 
around a fastener. This temperature can be 
subtracted from the total temperature change to 
calculate the thermal stress of a fastener with an 
initial clearance at room temperature. For a 
sample Joint the initial fastener fit is shown to 
have a significant effect on the thermal stresses. 

Using the algebraic equations developed in 
this study, material properties and Joint geometry 
were varied to determine their effect on thermal 
stresses. The effect of varying the radius of the 
structure surrounding the fastener was determined. 
When the ratio of structure to fastener radius 
exceeds five or six the peak stresses become 
insensitive to further Increases in the ratio. 

The effect of varying the structure-to-fastener 
modulus ratio was also shown. Increasing fastener 
modulus increases thermal stresses in the Joint. 
Thermal stresses of similar hollow and solid 
fasteners were compared. Hollowing out a fastener 
increases fastener stresses and decreases 
structural stresses. The effects of varying the 
Poisson's ratios of both the fastener and 
structure are also shown. Increasing the fastener 
Poisson's ratio effectively stiffens the fastener, 
which is in two-dimensional hydrostatic 
compression, and thereby increases the peak 
thermal stress. Increasing the structure 
Poisson's ratio, however, decreases peak thermal 
stresses. 

Finite element analysis was used to bound the 
effects of friction on thermal stresses and to 
show the effects of a fastener protruding beyond 
the surround structure. Friction is shown to have 
the potential to significantly Increase the 
thermal stresses in a Joint compared to the 
stresses calculated neglecting friction. The 
fastener protruding beyond the surrounding 
structure is shown to produce slightly higher 
localized stresses and to introduce relatively low 
out-of -plane stresses for a particular example. 

The equations developed in this paper are 
Intended to provide a designer with a quick. 
Initial estimate of the thermal stresses in and 
around a cylindrical fastener at elevated 
temperature. The equations should be used with a 
good understanding of the assumptions and 
limitations which went into their derivation. 


Appendix A 

Literature Review 

The problem of calculating thermal stresses 
around a cylindrical fastener can be reduced to an 
axlsymmetric, plane-stress problem if both the 
fastener and surrounding material are linearly 
elastic and isotropic. A general solution for 
such problems is presented by Little in reference 
2 and by Timoshenko and Goodier in reference 3. 


A similar problem is solved by Muskhellshvlli 
in reference A. In this problem an elastic 
circular disk is Inserted into a smaller radius 
hole in an infinite plate. This solution only 
gives the normal stress at the boundary between 
the materials and does not account for a finite 
annulus around the disk. Also, the solution does 
not involve the thermal expansion of the two 
materials. 

Several solutions for thermal stresses in the 
vicinity of cylindrical inclusions were found in 
the literature. References 5 and 6 present 
thermal stress sol utio ns fo r an infinite cylinder 
imbedded in an - elastic - sedium with different 
stiffness properties and thermal expansion 
coefficient in the presence of thermal gradients. 
The plane-strain solution for a uniform 
temperature change is presented in reference 5. 
When converted to a plane-stress solution the 
equations in reference 5 agree with the present 
solution for a cylinder with an infinite radius of 
material surrounding the cylinder. Similarly, the 
solution in reference 6, when simplified for a 
constant temperature change and converted to plane 
stress, also agrees with the current solution and 
that of reference 5 for a cylinder with an 
infinite radius of material surrounding the 
cylinder. The solutions in references 5 and 6, 
however, do not account for a finite radius of 
material surrounding the cylindrical inclusion. 

In references 7 and 8 a plane-stress solution for 
thermal stresses in a disk with a central shaft is 
given. The central shaft, however, is rigid and 
does not expand with temperature. A solution for 
thermal stresses in a circular plate with a 
central hot spot is also presented in reference 8. 
Reference 9 a more complicated problem with a pair 
of circular inclusions and a nonuniform 
temperature distribution. 

Although the general solution and similar 
solutions were found in the literature, no 
explicit solution for the thermal stresses around 
a cylindrical fastener was found. Consequently, 
algebraic equations for the stresses in a finite 
ring of material around a cylindrical fastener of 
a different material were developed in this study. 


Appendix B 

General Solution For Axlsymmetric Plane Stress 

The general plane-stress solution for 
axlsymmetric problems, as outlined in references 2 
and 3, is shown in the following mathematical 
development. The biharmonic equation which 
defines the Airy stress function for axlsymmetric 
plane stress is: 



Equation (B1 ) has a solution in terms of Integer 
powers of r: 

♦ - r n (B2) 

Substituting equation (B2) into equation (B1) 
produces a characteristic equation which leads to 
the general solution 
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♦ - A ln(r) ♦ B r* ln(r) ♦ C r ♦ D 


(B3) 


u - 


A (l*v) 


♦ 2 C (1-v) r 


In polar coordinates, 'the stresses in terns of the 
Airy stress function are 


o 


r 


1 3 ^ + _ 3^0 

r 3r r 1 2 36* 


♦ B (1-v) [2 r ln(r) -r]-2vBr+ const ) (BIO) 

u - g ( - - I*- - * 2 C (1-v) r 
p b r 

♦ B (1-v) [2 r ln(r) ♦ r] + 2 B r (B11) 


o 


6 




(BU) For equations (BIO) and (B11) to be consistent the 

constant B and the constant in equation (BIO) must 
be zero. The general expressions for plane stress 
are therefore: 


Therefore, by substituting the expression for 0 
into equations (B4), the stresses for an 
axisynmetric problem are found to be 

o « — 2 ♦ 2 B ln(r) ♦ B ♦ 2 C 
r r 

o Q - * 2 B ln(r) ♦ 3 B ♦ 2 C (B5) 

6 P 


0 


re 


0 


In polar coordinates, the definitions of strain 
are given by 



r 36 


(B6) 


° e - - Jr * 2 c 

and the radial displacement is given by: 

u r " f ( - * y— ♦ 2 C (1-v) r ) 

The constants A and C can be obtained by 
substitution of boundary conditions into 
(B12) and (B13). 


(B12) 


(B13) 


equations 
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